Stability of magnetic configurations in nanorings 
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The relative stability of the vortex, onion and ferromagnetic phases in nanorings is examined as 
a function of the ring geometry. Total energy calculations are carried out analytically, based on 
simple models for each configuration. Results are summarized by phase diagrams, which might be 
used as a guide to the production of rings with specific magnetic properties. 
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The understanding of the properties of small magnetic 
elements has been a major challenge in the rapidly evolv- 
ing field of nanoscale science. Besides the basic scientific 
interest in the magnetic properties of these systems, ev- 
idence is that they might be used in the production of 
magnetic devices, such as high density media for mag- 
netic recording Q, Q and magnetic logic 0, 0] • Common 
structures are arrays of nanowires cylinders 0, |(|, 
cones 0, and rings 0, 0. Recent theoretical stud- 
ies on such structures have been carried out aiming at 
determining the stable magnetized state as a function 
of their geometric details [ill Il2l PHil IT^j . Among the 
available geometries, the ring shape is of particular in- 
terest due to its core-free magnetic configuration leading 
to uniform switching fields, guaranteeing reproducibility 
in read- write processes 0. As a consequence, magnetic 
nanorings have received increasing attention over the last 
few years, from both experimental and theoretical points 
of view. 

Geometrically, ring-shaped particles are characterized 
by their external and internal radii, R and a, respectively, 
and height, H . Magnetic measurements as well as micro- 
magnetic simulations of such systems have identified two 
in-plane magnetic states, namely the flux-closure vortex 
state, V, and the so-called "onion" state, O. The lat- 
ter is accessible from saturation and is characterized by 
the presence of two opposite head to head walls 0, Il6| . 
In addition, for sufficiently high values of H, the occur- 
rence of ferromagnetic order, F, along the ring axis, is 
also possible. It is therefore clear that, for practical ap- 
plications, the determination of the ranges of values of 
the parameters R, a, and H, within which one of those 
configurations is of lowest energy, is of great relevance. 
However, magnetic measurements do not always allow a 
clear identification of the magnetic arrangement within 
the particles, which makes the theoretical approach to 
the problem highly desirable. In the present work we 
report results of total energy calculations for magnetic 
nanorings, on the basis of which the relative stability of 
the three above mentioned configurations could be deter- 
mined. Our results are summarized by phase diagrams 
in the R — H plane. 



Nanorings in the size range currently produced may 
consist of more than 10 8 magnetic atoms. As a conse- 
quence, the determination of the configuration of low- 
est energy based on the investigation of the behavior of 
individual magnetic moments becomes numerically pro- 
hibitive. In order to circumvent this problem, we re- 
sort to a simplified description of the system, in which 
the discrete distribution of magnetic moments is replaced 
by a continuous one, characterized by the magnetization 
M(r), such that M(f)SV gives the total magnetic mo- 
ment within the element of volume 5V centered at r. We 
also require that |M(f)| = Mq, the saturation magneti- 
zation. 

The internal energy, E to t, of a single ring is given by 
the sum of three terms corresponding to the magneto- 
static (Edip), the exchange (E ex ), and the anisotropy 
(Ek) contributions. Since nanorings are usually poly- 
crystalline, the magnetic anisotropy averages to zero due 
to the random orientation of the crystallites |17| . In view 
of that, it will be neglected in our calculations. 

The dipolar contribution can be obtained from the 
magnetization according to the relation 



Edip — 



VU(f) dv 



(1) 



where an additional configuration independent term has 
been left out, and the magnetostatic potential is given 
(in SI units) by 



r — r 47r 



M{f 



■ds'. (2) 



In the above expression, V and S represent the volume 
and surface of the ring, respectively. Assuming that the 
magnetization varies slowly on the scale of the lattice 
parameter, the exchange energy can be written as 



E r 



A 



(Vtoz) 2 + (Vm y ) 2 + (Vm,)' 



dv 



(3) 
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where mj = Mj/Mo (i = x,y,z) are the reduced com- 
ponents of the magnetization and A is the stiffness con- 
stant of the material, which depends on the exchange 
interaction between the magnetic moments. Also in 
this expression, an additional configuration independent 
term has been left out 0|- At this point it is conve- 
nient to define dimensionless variables. The natural scale 
for linear dimensions is given by the exchange length, 
L ex = y // 2A/ hoMq , in terms of which we can define 
the dimensionless radius and height p — R/L ex and 
h = H/L ex , respectively, ft is also worth defining the 
parameter 7 = H/R and the aspect ratio [3 — a/R, such 
that < /3 < 1. 

In order to proceed with our calculations, it is neces- 
sary to specify the function M(r) corresponding to each 
of the three configurations under consideration. For the 
ferromagnetic state, M(f) can be approximated by M z, 
where z is the unit vector parallel to the axis of the ring, 
whereas for the vortex state we can take M(r) = Mo</>, 
where 4> is the azimuthal unit vector in the ring plane, xy. 
The most interesting case, however, regards the onion 
configuration, whose typical arrangement of the magnetic 
moments is illustrated in Fig.(l). So far, no analytical 
model has been presented in the literature to describe 
such configuration. However in a recent work, Beleggia 
et al. lj| have investigated the magnetic phase diagram 
for rings, considering a ferromagnetic-in-plane configura- 
tion instead of the onion configuration. 




any model for the onion configuration has to allow the 
magnetization either to approach that of an in-plane fer- 
romagnetic or to deviate from it, depending on the values 
of these parameters. This point will be made clearer fur- 
ther on in this article. 

Taking into consideration the spacial symmetry of M, 
as represented in Fig.(l), the expressions for the reduced 
components m r — M r /Mo and m^> = M^/Mq of the 
magnetization can be written as 



m r ((f>) 



f ((/>), < </>< tt/2 

-f(n - 0), 7T/2 < <j) < 7T 

-f((f) -7T), 7T<0<37T/2 

/(2tt -</»), 37r/2<0<27r 



(5) 



and 



m,p((j)) 



- Vl-PW, O<0<^/2 
- yi-p(n-<t >), tt/2<^<tt 

VI - P{4> - g, 7T<0<37T/2 

^1-/2(2^-0), 3tt/2 < <f> < 2tt 



(6) 



where / is a bounded function such that — 1 < /(</>) < 1. 
A suitable form for the function / is 



/(n,0=cos n (0) , 



(7) 



where n is a continuous variable, defined such that n > 1, 
chosen so as to minimize the total energy of the config- 
uration for each value of p, h and (3. Such model allows 
a continuous transition from the in-plane ferromagnetic 
state (n = 1) to the onion state (n > 1). The deviation 
from the ferromagnetic configuration becomes more pro- 
nounced as n increases from 1. Fig.(l) shows a possible 
onion configuration with n > 1. 

Having specified the functional form of M(r) for the 
three configurations, we are now in position to evaluate 
the total energy for each of them. 



A. Ferromagnetic configuration (F) 



FIG. 1: Onion configuration of a Co ring with h = 2.1, p = 7 
and P = 0.5 for n = 4.6. The arrows denote the direction of 
the magnetization. 

The model we adopt is based on the assumption that 
the magnetization in the onion configuration can be ap- 
proximated by 



M(f) = M r {<P)r + M,p 



(4) 



where the possible dependence of M r and of on r has 
been neglected. The important point here is the implicit 
dependence of the two components of the magnetization 
on the geometry of the ring, more precisely, on the val- 
ues of p, h and (3. Therefore, in order to be realistic, 



From Eq.lO we immediately find that E^ x — 0. Thus, 
the total energy has only the dipolar contribution, which 
can be obtained from the magnetostatic potential U(r) 
given by Eq.Q. Using the expansion l|17|) in expression 
@ we obtain 



El 



MoM 2 i| 



P 3 *(/3,7) , 



where the function 7) is defined by 



*(/3, 7 ) 



(8) 



1 e 2 ^ [Ji (V) - P-h 0%)] 2 dy . (9) 
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Here Ji(z) are Bessel functions of the first kind. 

Equation JSJ has been previously obtained by Beleggia 
et al. (See Eq.(15), Ref. considering a more general 
approach based on Fourier Transforms of the magnetiza- 
tion. 



W,l) = - I dx I 6y[J p (y)-pJ P (py)] 



o 



e lv + jy — 1 

y 2 



J P (xy) , 



B. Vortex configuration (V) 

It is clear from Eq.(0l that for this configuration 
E^ ip = 0. The exchange contribution can also be evalu- 
ated from Eq.lO so that the final result for the reduced 
total energy reads 



E. 



E 



IM)MqL\ x 



-irh In [3 



(10) 



C. Onion configuration (O) 

In the previous two configurations, the magnetization 
function satisfies the condition V • M = 0, which means 
absence of volumetric magnetic charges. In the onion 
configuration, however, due to the presence of two regions 
with head-to-head domain walls, in which V - M ^ 0, 
the dipolar energy turns out to be given by the sum of 
two contributions, E® s and E® v , coming from the sur- 
face and volume terms in the expression for U(r), respec- 
tively. Details of the calculations of these two quantities 
are given in the Appendix. Results for the reduced ener- 
gies read 



E. 



dS 



P 3 J2 b p( b pQi + d pQ^ > ( n ) 

P =i 



h dV = 



E. 



dV 



P 3 y^( d p - b p) (hQs + dpQi) , 



P =i 



where 



tt/2 



b p (n) 



f(n, <j>) cos(p<f>) Acj) , 



(12) 



(13) 



7T/2 



d p (n) = p J Vl-/ 2 (n,0)sin(p^ deb , (14) 
o 



and 



?i(/?,7) = - I xdx I dy [J p {y) - pj p (py)} 



e -lV + 7 y - 1 
V 



[Jp-i(xy) - J p +i(xy)] , 



3 (A 7) = - / dw J x dx I dy [J p _i(xy) - J p+1 (xy)} 
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e -lV +jy-l 
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1 1 00 

Q 4 (/3,7) = ^y d w J J dy J p (wy) 



' e~ lv + jy — 1 
y 2 



J P (wy) , 



J p (xy) . 



We remark that the sum in Eqs. and (|12[) runs over 
odd values of p and converges quite rapidly so that, in 
practice, just a few values of p (say, 5 or 6) have to be 
considered. 

The exchange energy can be also obtained analytically, 
as indicated in the Appendix. The reduced energy Ef x 
is then given by 



F,° 



where 



I{n) 



tt/2 



-irh [I(n) - 1] In /? , (15) 



df(n, 0) 



1 - / 2 (n, 



As already pointed out, the value of n in the above 
expressions is chosen so as to minimize the total energy, 
for fixed values of the geometrical parameters (3, 7 and 
p. In other words, it is obtained by solving the equation 
dE°Jdn = 0, where E? ot = Eg+Ef y +Eg. It is worth 
pointing out that the dependence of Ef ot on n enters just 
via the functions b p (n), d p (n) and I(n). The first one can 
be expressed in terms of Gamma functions, whereas the 
remaining two can be easily evaluated numerically. On 
the other hand, the functions Qj can be calculated just 
once for given values of (3 and 7. 



III. RESULTS AND DISCUSSION 

The above expressions for the total energy enable us 
to investigate the relative stability of the three configu- 
rations of interest. For each aspect ratio (3 we can de- 
termine the ranges of values of the dimensionless radius 



p and height h within which one of the three configura- 
tions is of lowest energy. The boundary line between any 
two configurations can be obtained by equating the ex- 
pressions for the corresponding total energies. Figure (2) 
illustrates phase diagrams for different values of (3. The 
stability regions corresponding to each configuration are 
indicated in the diagram. 
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FIG. 2: Magnetic phase diagrams for rings with (3 = a/R — 
0.2 (diamonds), 0.5 (full circles), and 0.8 (stars). F, V and O 
indicate regions where ferromagnetic out-of-plane, vortex and 
onion configurations are relatively more stable, h and p are 
the dimensionless geometric parameters, defined respect to 
the material exchange length as h — H/L ex and p — R/L ex . 
The inset illustrates the geometry of a ring. 

It is interesting to look at the dependence of the expo- 
nent n, which determines the magnetization profile in the 
lower energy onion configuration, on the geometry of the 
ring. Fig. (3) presents a set of lines in the p/i-plane corre- 
sponding to constant values of n, for (3 — 0.5. The points 
on such lines laying outside the O region (open symbols) 
correspond to cither metastable or unstable onion con- 
figurations. We shall come back to this point further on 
in this article. We clearly see that for small rings, the 
onion configuration turns out to be rather close to the 
in-plane ferromagnetic one [n = 1). Another interesting 
point is the weak dependence of the in-plane magnetiza- 
tion of the onion configuration as a function of n. The 
reduced magnetization rn x can be obtained by integrat- 
ing M x within the volume of the ring, which gives us 
m x = (2/7r) [bi(n) + d 1 (n)]. The inset in Fig. (3) presents 
function of n, which well illustrates this point. 

The transition between the F and O configurations is 
determined by the balance between the energies of the 
out-of-plane and of the in-plane magnetic configurations. 
A qualitative understanding of such transition can be 
achieved by examining the strength of the demagneti- 
zation fields along the z and x axes. For given (3 and 
p, larger values of h result in reduced demagnetization 
fields along the z direction, favoring the F configuration. 
This is the case, for example, of rings with p = 1.0 and 




FIG. 3: Constant n curves for rings with j3 = 0.5. Solid 
symbols correspond to stable onion configurations, and open 
symbols to metastable or unstable ones. Thick solid lines 
denote the boundaries between stability regions and thin lines 
are guides to the eye. Results are presented for n = 1.01 
(triangles), n = 1.1 (squares), 1.3 (circles), 1.5 (diamonds), 2 
(stars), and 4 (pentagons). Inset: magnetization component 
along the onion axis (m x ) as a function of constant-n curves in 
the p/i-plane for rings with (3 = 0.5. Results are presented for 
rings both inside (full symbols) and outside (open symbols) 
the onion stability region. The latter correspond to either 
unstable or metastable onion configurations. 



(3 = 0.2, when h is increased from 1.0 to 3.0. On the 
other hand, for fixed p and h, a decrease in [3 is equiva- 
lent to a reduction of the inner radius, which results in 
a larger in plane demagnetizing field. Thus, a ring with 
p = 1.0 and h = 1.0 exhibits an in-plane magnetic order 
for (3 = 0.2, and an out-of-plane order for (3 = 0.8. 

Concerning the two in-plane configurations, namely 
O and V, the transition between them depends on the 
strength of the demagnetization field along the x direc- 
tion in the domain wall regions of the O configuration. 
The larger the value of the ring width W = R — a, the 
smaller this field is, favoring the onion state. Thus, a 
ring with p = 3.0 and h = 1.0 exhibits the V phase for 
(3 = 0.8, and the O configuration for (3 — 0.2. 

The mechanism responsible for the transition between 
the V and F configurations is somewhat subtle. We no- 
tice in Fig. (2) that the line separating these two regions 
exhibits a maximum shift to the left for [3 about 0.5. 
As a consequence, rings with say p = 3.0, h = 6.0 and 
inner radius corresponding to (3 = 0.2, 0.5, and 0.8 are 
found to exhibit F, V , and F configurations, respectively. 
Such interesting behavior can be understood by consid- 
ering the relative strength of the exchange and dipolar 
energies in the two states. It is clear from Ea. (|10|) that 
for fixed values of p and h the exchange energy in the 
vortex state diverges at (3 — 0, whereas that for the F 
configuration approaches a finite value [cf. Eqs© and 
JUJ], Thus, for sufficiently small inner radii (f3 <C 1), the 
F configuration has lower energy than the V one. The 



5 



reason for that is the large contribution to E^ x coming 
from the central region of the ring. Then, the increase in 
the inner radius leads to a rapid decrease of E^ x , reducing 
EY ot with respect to Ef ot and favoring the V configura- 
tion for intermediate values of (3. However, for values of 
(3 close to 1 (narrow rings) and sufficiently large values 
of h, the dipolar energy in the F state, which can be as- 
sociated to the magnetic charges at the top and bottom 
surfaces of the ring, becomes smaller than E^ xl favoring 
the ferromagnetic state. 
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FIG. 4: Triple point position (p t , ht) as a function of j3 = a/R 
(full circles). The result for a cylinder is represented by a star. 

All diagrams in Fig. (2) show a triple point (p t ,ht), 
corresponding to the situation in which the three config- 
urations have equal energy. Fig. (4) shows the trajectory 
of the triple point (r t (/3), h t (/?)) as a function of (3. We 
notice that the curve correctly extrapolates to the result 
for cylinders, which has been numerically obtained by 
d'Albuquerque e Castro et al. [2(j using a scaling tech- 
nique, and analytically, by Landeros et al. [2l|. 

An interesting result is the occurrence in some cases of 
metastable onion configurations. Castaho et al. |l6j ] have 
carried out a detailed experimental study of the magnetic 
behavior of Co nanorings with H = 10 nm and external 
and inner radii ranging from 90 nm to 260 nm and from 
20 nm to 140 nm, respectively. For rings with R = 180 
nm and a = 20 nm, and R = 180 nm and a = 70 nm, 
they have found that the corresponding hysteresis curves 
clearly indicate the occurrence of metastable onion con- 
figurations. Such states can be reached by reducing the 
applied magnetic field from the saturation value to zero. 
However, as the magnetic field is increased in the op- 
posite direction and the magnetic energy becomes suf- 
ficiently high to overcome the energy barrier between 
the O and V configurations, the systems undergoes a 
transition to the latter state. Such behavior is entirely 
consistent with our results. Indeed, taking into consid- 
eration that for Co L ex — 2.85 nm, the rings examined 
by Castaho and collaborators have reduced dimensions 
p = 63.2, and h = 3.5 and correspond to (3 = 0.11 and 
0.39. Thus, from the results in Fig. (2), we expect the 
two rings to fall well within the vortex region. We re- 



mark that from Fig. (4), the upper bound for the posi- 
tion of the triple point for f3 — > 0, is p max = 3.8 which 
is much smaller than the value corresponding to the two 
samples. As a consequence, the observed onion configu- 
rations in such systems represent necessarily metastable 
states. The occurrence of metastable onion states within 
the vortex region has also been found in micromagnetic 
calculations. Using the OOMMF package [2^ we have 
carried out ground state calculations for a Co ring with 
h = 2.1, p = 7, and (3 — 0.5. We have found that the 
system might be trapped in a metastable onion configu- 
ration when the in-plane ferromagnetic state is taken as 
the starting configuration. These results are presented in 
Fig.(l). 

In conclusion, we have theoretically investigated the 
dependence of the internal magnetic configuration of 
nanorings on their geometry. Three typical configura- 
tions have been considered, namely out-of-plane ferro- 
magnetic, vortex, and the so-called onion configuration. 
For the latter, we have proposed a simple analytical 
model, which allows a continuum transition between the 
onion and the in-plane ferromagnetic states. Our results 
are summarized in phase diagrams giving the relative sta- 
bility of the three configurations. The possibility of the 
systems assuming unstable or metastable configurations 
has also been considered. Our results might be used as a 
guide to experimentalists interested in producing samples 
with specific magnetic properties. 
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IV. APPENDIX 

The calculation of the dipolar energy of the onion con- 
figuration begins by replacing the functional form Eq. Q 
in Eq. leading to 



Mo 
2 



M r (<j>) 



dU + M {cj)) dU 
dr r dd> 



dv . 



(16) 



6 



For the calculation of the magnetostatic potential we use Using 
expression |QJ and the following expansion [23| 



oo 00 

_L_ = J- e*W-*'> J J p (kr)Jp(kr')e k ^<'^dk 1 



P——00 



where J p are the first kind Bessel functions. This way 
the surface contribution to the potential (second term in 
Eq.@), reads 



TI Mo V 



2tt 



P——00 



H 



dfc Jp(kr) [RJ p (kR) - aJ p (ka)} / e k{z <- z > ] dz' 



Using m r ((f) defined by Eq. JSJ it is straightforward to 
obtain 



271 



?7i 



r (cf) e"** d0 = 2 (1 - e"**) 6 p , 



where b p is defined by Eq. This relation lead us to 

consider only odd values of the sum index p in what fol- 
lows. Integrating over z', and after some manipulations, 
we obtain that 



00 



Us — ^ b p cos(p(f>) I dfc J p (kr) 

2 — e~ kz — e - k i H ~ z 



x [i?J p (fci?) — aJp(ka) 
Introducing this potential in Eq. (|16f) we obtain 

00 H 

^ = / ^ &/ d^ 



7T ^ P J k 

P =1 



2tt 



m r (</>) cos(p(f) d0 = ib p 



(17) and 



2/r 



p m<j, ((f)) sin(p0) d</> = -4d p , 



in Eq. 1|18[1 we obtain 
£ d ° 5 = ^MlR^l J d2/ 



P =i 



[J P (y) - PJp(Pv)] 



e -iv +7^-1 
2/ 2 



da; |6pd p J p (xy) + 6 p ^ [J p _i(xy) - J p+i (a;y)]| 



Thus, the superficial contribution to the reduced dipolar 
energy (Eq. JTIJ) can be written as 



0=1 



2 — e ^ fc ^ — e -k(H-z) 



The calculation of the volumetric contribution (Eq. 1)12(11 
follow the same procedure. 

The exchange energy of the onion configuration comes 
from substituting the Cartesian magnetization compo- 
nents 

m x (4>) = m r ((f>) cos((f) - m0 ((f)) sin(</>) 
m y ((f) = to,, ((f)) sin(</3) + ((f) cos(</>) , 

in the semiclassical expression Eq. from which we 
obtain that 



R 



kr 



x/ dr [RJ p (kR) - aJ p (ka)]{ — [J p _i(fcr) - J p+1 (kr)} 



xj d(f [m r ((f) cos(p(f) — Jp(kr)pm t j } ((f)) sm(p(f)] } . (18) 




^,r> h , 1 



fdm r Y fdm^ 



By using Eqs. and JBJ) we found the reduced exchange 
energy of the onion configuration (Eq. I15fl . 
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